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An experimental study on the design of counter-rotating axial-flow fans was carried out. The fans 
were designed using an inverse method. In particular, the system is designed to have a pure axial 
discharge flow. The counter-rotating fans operate in a ducted-flow configuration and the overall 
performances are measured in a normalized test bench. The rotation rate of each fan is independently 
controlled. The relative axial spacing between fans can vary from 17% to 310%. The results show 
that the efficiency is strongly increased compared to a conventional rotor or to a rotor-stator stage. 
The effects of varying the rotation rates ratio on the overall performances are studied and show 
that the system has a very flexible use, with a large patch of high efficient operating points in the 
parameter space. The increase of axial spacing causes only a small decrease of the efficiency 



I. INTRODUCTION 

Early studied in the 1930's 1 -, the counter-rotating machines arouse a greater interest in the turbomachinery 
field, particularly for their potential improvement of the efficiency with respect to conventional machines 
by recovering kinetic energy from the front rotor exit-flow and by adding energy to the flow. The first 
counter-rotating machines have appeared in aeronautic 1 -^ and marine applications 5,4 in open configuration. 

Conventional designs of high speed counter-rotating fans are based on quite expensive methods and 
require a systematic coming and going between theoretical methods - such as the lifting line theory or 
the strip-analysis approach^- and CFD analysis^. Moreover, the axial spacing, which has a major role on 
the rotors interaction and consequently on the noised, is a key parameter to find a compromise between 
high aerodynamic and good acoustic performance for high speed fans 5 . In order to reduce this interaction, 
the axial spacing of high speed fans has to be relatively large, resulting in a decrease in the aerodynamic 
performance 5 -. For the same reason, the rear rotor (RR) diameter has to be smaller (about 10% according 
to 5 -^) than the front rotor (FR) diameter to reduce interaction between the FR tip vortex and the RR blade 
tip. 

Contrary to that, in the case of low speed fans axial spacing could be shortened using the benefit of a 
relatively low rotor interaction. Therefore these machines see a revival of interest in several distinct configu- 
rations -open and ducted flows, shrouded or not shrouded rotors- in various subsonic regime applications^—. 

Recent research work dealt with the effects of global parameters like rotation speed ratio 1 ^, local phenom- 
ena such as tip vortex flowa 1 ^ and improvement of cavitation performance for pumps^. All previous studies 
have shown the benefit of RR in improving the global efficiency and in increasing the operating flow-rate 
range while maintaining high efficiency. The counter-rotating systems (CRS) moreover allow to reduce the 
fans diameter and/or to reduce the rotation rate. More axial spacing is needed compared to one simple fan, 
but not much more than a rotor-stator stage. However, it requires a more complex shaft system. Another 
interesting feature of CRS is that it makes it possible to design axial-flow fans with very low angular specific 
speed = ( A ^/^3/4 with oj = ""'° rl +" r °'° r:i the mean angular velocity, Q the flow rate, Ap t the total 
pressure rise, and p the fluid density. With such advantages, the CRS becomes a very interesting solution 
and the interaction between the rotors needs to be better understood in order to design highly efficient CRS. 
However, only a few studies have been concerned with, on the one hand, the effect of the axial spacing, and, 
on the other hand, the design method-^, particularly with rotors load distribution for a specified design 
point. 

This paper focuses on two major parameters of ducted counter-rotating axial-flow fans in subsonic regime: 
the rotation rate ratio, 6 and the relative axial spacing, A. In some cases, these systems are studied by using 
two identical rotors or the RR is not specifically designed to operate with the FR. In this study, the FR is 
designed as conventional rotor and the RR is designed on purpose to work with the FR at very small axial 
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spacing. In this first design, The total work to perform by the CRS was arbitrarily set up approximately to 
two halves one half respectively for the FR and RR. In § |TT] the method that has been used to design the 
front and the rear rotors is firstly described. The experimental set-up is presented in § IIIII Then the overall 
performances of the system in its default configuration and the effects of varying the rotation ratio and the 
relative axial spacing between the rotors are discussed in § IIVI 



II. DESIGN OF THE ROTORS 



A. General approach 



The design of the rotors is based on the use of the software MFT (Mixed Flow Turbomachinery) , a ID code 
developed by the DynFluid Laboratory^ — based on the inverse method with simplified radial equilibrium — 
to which an original method has been added specifically for the design of the RR of the counter-rotating 
system. 

From the specified total pressure rise, volume flow-rate and rotating speed, optimal values of the radii 
Rup and Rhub are first proposed. In a second step, the tip and the hub radii as well as the radial distribution 
of the circumferential component of the velocity at the rotor outlet, C U 2(r), could be changed by the user. 
The available vortex models are the free vortex (C U 2(r) = ^r), the constant vortex (C„2(r) = K ) and the 
forced vortex (C U 2(r) — rK ). 

The velocity triangles are then computed for 11 radial sections, based on the Euler equation for perfect fluid 
with a rough estimate of the efficiency of rj es t = 60% and on the equation of simplified radial equilibrium 
(radial momentum conservation). The blades can then be defined by the local resolution of an inverse 
problem considering a 2D flow and searching for the best suited cascade to the proposed velocity triangles 
by the following parameters: 7 the stagger angle, computed from the incidence angle, a giving the lower 
pressure variation on the suction surface of the blade using equations [1] and [5] The solidity, a and the chord 
length, c are thus computed at the hub and at the tip using equations [4] and [5] where V denotes the Lieblein's 
diffusion factor—. The intermediate chords are obtained by linearisation. Finally, the camber coefficients 
CzooO are computed using equation [51 

7 = ft -a (1) 

g(ft) = 2.103 -4.01910" 7 ^f 382 (3) 
*-l + £W^ (4) 
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a + 2.525 

C zoo0 = r~, 0.823 6 

p{<j) = 15.535 - 12.467e- 4242CT (7) 

These empirical equations have been validated for NACA-65 cascades^, for 0.5 < a < 1.5 and < C zoo o < 
2.7. 

The behaviour of the designed machine resulting from the above method can then be analysed using a 
direct method in order to determine whether the design point is achieved and what are the characteristics 
of the machine at the neighbourhood of the design point. The effects due to real fluid are taken partially 
into account with in-house loss models and the introduction of an axial- velocity distribution which considers 
the boundary layers at the hub and casing. Thus, the characteristics of the machine can be obtained in the 
vicinity of the design-point discharge. 

Regarding the CRS, the geometrical dimensions, the number of blades of FR and of RR and their rotation 
rates are imposed. In particular, the number of blades of each rotor was chosen in order to prevent to have 
the same blade passing frequency or harmonics for both rotors in the lower frequencies range. The system 
that is presented here has moreover been designed to have a pure axial exit-flow. An iterative procedure is 
then performed. The pressure rise of the FR is initially chosen and then designed and quickly analysed as 
explained. An estimate of the pressure rise that RR would made is then performed, based on this analysis. 



FIG. 1. Velocity Triangles for the CRS. The fluid is flowing from left to right. 



TABLE I. Design point of the counter-rotating system for air at p — 1.21 kg.m 





CRS 


FR 


RR 


D (mm) 


380 


380 


380 


Ru P (mm) 


187.5 


187.5 


187.5 


Rhub 

(mm) 


55 


55 


55 


Z 




11 


7 


Apt (Pa) 


420 


260 


160 


TV (RPM) 


1900 


2000 


1800 


Q 

(m 3 .!!" 1 ) 


3600 


3600 




n 


2.46 


3.71 




Other con- 
straints 


Axial exit- 
flow 


Constant 
vortex 





If the total pressure rise of the CRS is not met, the design pressure rise of FR is varied and the calculus are 
made again. In this method, losses and interactions in-between the two rotors are not taken into account. 
Any recirculation happening near the blade passage or near the blade hub or tip is not predicted by MFT 
as it is based on simplified radial equilibrium. 



B. Design of the Front Rotor 

The design point of the CRS is given in Tab. Q] The system is designed to achieve a total-pressure rise 
Ap t = 420 Pa at flow-rate Q — 3600 m 3 .h. _1 for a rotation rate around 2000 rpm (ft ~ 2.46). The geometrical 
dimensions are fixed to fit in the set-up: Ru p — 187.5 ram and *j* ub = 0.293. In this first and simple design, 
the front rotor (FR) has arbitrarily larger total pressure rise than the RR. The constant vortex model leads 
to a uniform tangential velocity distribution downstream FR for perfect fluid and thus has been used to 
simplify the design of the RR. Also, the rotors are not shrouded and the radial gap between the blade tip 
and the wall casing is of 2.5 mm, i.e. 1.9% of the blade height. 

In the present case, the direct analysis predicts a mean absolute tangential velocity at the design flow-rate 
Cq2FR — 9.6 m.s -1 with a radial distribution uniform within ±5% (constant vortex design). The Reynolds 
number based on the inlet relative velocity and the chord varies from 0.6 x 10 5 at the hub to 3 x 10 5 at 
mid-span and 7 x 10 5 at the tip. 
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TABLE II. Blade cascade parameters for the two rotors. Radius R (mm). Chord length c (mm). Cascade solidity 
a. Stagger angle 7 (°). Profile designation according to the nomenclature given in Ref.J**: NACA65(xx)yy with (xx) 
representing the relative camber and yy standing for the relative thickness. Lieblein's diffusion factor T> 



R c a 7 profile T> 


Front Rotor (blade thickness 4.5 mm) 


Hub 


55 40.31.2823NACA 65(26)11 0.62 


Mid- 


121. 2558.00. 8457NACA 65(12)07 


span 




Tip 


187.5 75.70.7169NACA 65(07)06 0.44 


Rear Rotor (blade thickness 6 mm) 


Hub 


55 58.81. 1 873 N AC A 65(03)10 0.61 


Mid- 


121.2572.90.6665NACA 65(05)08 


span 




Tip 


187.5 87.10.5175NACA 65(04)07 0.46 



C. Design of the Rear Rotor 

The method used for the design of the RR is to consider the velocity and the flow angle at the trailing 
edge of the FR blades. Therefore, FR was analysed with MFT to retrieve the axial and tangential velocities 
(C a 2FR — C a iRR an d C U 2fr. = V u iRR respectively) and the angle ct2R\ in the absolute reference frame, at 
the exit and along the blade as shown in Fig. [TJ Therefore, the Euler work distribution along the blade does 
not match with any of the vortex models, previously mentioned. 

Using the same radial inverse design equations on 11 radial sections, the RR is drawn in such a way that 
the exit flow is purely axial, that is a.2RR(r) = 0°. The second hypothesis is that the axial velocity profile 
is kept constant across RR, i.e. C a 2RR{r) = C0 a \RR(r). Under these assumptions, the total pressure rise 
of RR should be Ap t RR — r/ est pU m RRC U 2FR — 0.6 x 1.2 x 22.9 x 9.6 ~ 160 Pa where rj es t = 0.6 is an 
empirical estimated efficiency observed from previous experimental designs. The blade cascade that lead to 
the desired velocity triangles is then designed with the previously described inverse method, adjusting the 
free parameters in such a way that the solidity lays in the range 0.5 < a < 1.5 and that the camber lays 
in the range < C 200 o < 2.7. After several iterations, the RR was drawn with Z = 7, T>h u b = 0.61 and 
Du p = 0.46. The geometrical characteristics of the rotor blades obtained with this method are summarized 
in Tab. HU Pictures of the Front and Rear rotors are given in Fig. [3J 



III. EXPERIMENTAL SET-UP 
A. Test bench 

The counter-rotating system is studied in a ducted-flow test rig — AER0 2 FANS— that has been built 
according to the ISO-5801 standards (installation of category B) 17 ' 18 . The Figure [5] shows this test rig. Two 
brushless PANASONIC A4 motors drive each rotor separately and are hidden in a casing of diameter 0.33D 
and of length 0.45D, with a warhead-shape end. The front and the rear motors are bound to the tube by 
two rod rows (3 and 5 rods, the first row being at 0.1D from the RR). For the front motor a honeycomb 
is placed between the two rows to reduce the rods effect on the incoming flow at the inlet of the FR and 
to homogenize the inlet flow. An iris damper — originally used for air flow regulation in ducts — is placed 
at the exit of the pipe to vary the test-bench hydraulic impedance and thereby to vary the operating point 
of the studied axial-flow fan. Finally, an axial blower can also be used at the exit of the pipe to widen the 
explored flow-rate. 




FIG. 2. Top: picture of the front rotor (left) and Rear Rotor (right). Bottom: experimental facility for CRS, 
AERO 2 FANS 
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TABLE III. Nominal points of FR rotating alone at Nfr = 2000 rpm, RR rotating alone at Nrr — 1800 rpm and 
CRS at N FR = 2000 rpm and 9 = 0.9 (see also Fig. 





FR 


RR 


CRS 


Max efficiency 46.2±1% 
(%) 


54.4±1% 


65.1±1% 


Nominal 
(m 3 .!!- 1 ) 


Q 3636±36 


2520±36 


3600±36 


Ap s (Pa) 


157±3 


88±3 


335±5 



B. Measurements methodology 



The study focuses on the influence of the relative axial spacing A — where S and CpR m are the axial 

spacing and the FR chord length at mid-span respectively, as well as on the influence of the rotation rate 
ratio 9 = ^ RK - Six axial spacings, from A = 0.17 to A — 3.1 are investigated by means of blocks of different 
thickness, "file reason of positioning RR the closest to FR and then increasing A is to investigate on any 
possible potential effects due to the rotors proximity. Regarding the rotation rate ratio, each rotor is driven 
separately so all combinations are possible and the only limitation is the motor power. Unless specified the 
default axial spacing is A — 0.17 and the default rotation rate ratio is 9 — 0.9 (see Tab. [XJ) . 

The static pressure rise of the rotor is obtained according to the ISO-5801 standards by measuring the 
pressure difference between the atmosphere and downstream the anti-gyration device. The static pressure 
losses mainly induced by the honeycomb, the motors casings and the anti-gyration device have been measured 
using an auxiliary axial blower (with both rotors removed) and have been added to the static pressure rise 
measurements . 

In the design of the CRS, it has been imposed a pure axial- flow at the exit of RR. In that case the static 
pressure rise of the CRS should be Ap s = Ap t - 1/2 p (Q / (irD 2 /4)) 2 ~ 373 Pa. If the flow at the exit is 
not purely axial, then the total pressure rise will remain the same but the static pressure rise will be smaller 

-. 2 

and equal to: Ap s ~Ap t — ( 1/2 joC^r) ). 
The static efficiency is defined by equation [5] 

^ = AP S Q (8) 

{Tfrufr) + (TrrOJrr) 

The torque T was measured using the drivers provided with the motors. A calibration measurement has 
been performed with a conventional torque-meter. This calibration step shows that the torque supplied by 
the driver is very close to that given by the torque-meter (relative error of 0.5%). 

Finally, for all performance measurements, the atmospheric pressure, the dry temperature and the wet 
temperature were measured and thus the density was computed for each measurement. It has been found 
that the variation in density, relatively to the design density, p air = 1.21 kg.m~ 3 is between 0.5% and 2.2%. 
Therefore, in order to present homogeneous results, the pressure rise is multiplied by the ratio of design 
density p a i r over the experimentally measured density, p exp , i.e. £ait. 

Pexp 



IV. RESULTS AND DISCUSSION 



A. Overall performances of the reference system ({8 = 0.9; A — 0.17}) 

The characteristics of the FR rotating alone (RR has been removed from its shaft in that case), of the 
RR rotating alone (FR has been removed) and of the counter-rotating system are shown in Fig. |3l 

The nominal flow-rates of the three systems, i.e. the flow-rates at maximum efficiency, are reported in 
Tab. Mil together with the corresponding static pressure rises and efficiencies. 

The FR rotating alone has a very flat curve (□ in Fig. [3]). The nominal flow-rate of FR is slightly greater 
than the design point — it is 3% greater. The measured static pressure rise at the design point is 157 ±3 Pa, 
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FIG. 3. Fans characteristics: (a) static pressure rise Ap s vs flow rate Q; (b) static efficiency 7/ s vs flow rate Q. The 
axial spacing is A = 0.17. □: FR rotating alone at Nfr — 2000 rpm (RR has been removed), RR rotating alone 
at N RR = 1800 rpm (FR has been removed) and o: CRS at N FR = 2000 rpm and 6 = 0.9. The ■ and the dashed 
lines stand for the design point of the CRS 



with a relatively low static efficiency of 46.2%. This is not surprising with no shroud and a large radial gap. 
Moreover, this is consistent with the estimated static pressure rise by MFT, which is around 151 Pa. 

The RR rotating alone has a steeper curve in Fig. [3]) and its nominal flow-rate Q ~ 2600 m 3 .h _1 
is lower than the design flow-rate of FR and CRS. This is consistent with the bigger stagger angle of the 
blades (see Tab. HI)) and can be explained by examining the velocity triangles in Fig. [T] and considering the 
case with the FR coupled to the RR: the incoming velocity C\rr=Cifr has an axial component as well as 
a tangential component. Hence, the flow angle in the relative reference frame reads: 



tan(/3 1RR ) 



URB, + CulRR 



Col 



RR 



(9) 



Now the case without the FR is considered and it is assumed that the flow through the honeycomb is axial. 
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Since the tangential component does not exist any more, C u \rr — m.s^ 1 . MFT estimates (Urr) ~ 
22.9 m.s" 1 , {Caiim) — 8.8 m.s -1 and (C u \rr) = {C u2 fr) — 9.6 m.s -1 , which leads to (tan((3inu)) ~ 3.69 
at the blade mid-span. Supposing that RR rotating alone reaches its maximum efficiency for (tan(/3i RR )) ~ 

3.69, equation [9] implies that (C a i RR ) = t an^(pf RR )) - 6 - 2 m - s_1 ) Le - Q - 2540 m 3 -h _1 . This is exactly the 
nominal flow-rate of RR rotating alone (see Fig. [3] and Tab. IIII[) . It is clear from the above analysis why the 
nominal flow-rate of RR is lower than the design flow-rate. 

The characteristic curve of the CRS (o in Fig. [3]) is steeper than the characteristic curve of FR. It is 
roughly parallel to the RR curve. The nominal flow-rate of the CRS matches well with the design flow-rate, 
i.e. 1 m 3 .s _1 . The static pressure rise at the nominal discharge (Ap s cns — 335 Pa) is 10.2% lower than the 
design point (373 Pa), which is not so bad in view of the rough approximations used to design the system. 
Please notice that the static pressure rise of the CRS is not equal to the addition of the static pressure rise 
of the FR with the pressure static rise of the RR, taken separately. 

The CRS has a high static efficiency {r) s CRS — 65%) compared to a conventional axial-flow fan or to a 
rotor-stator stage with similar dimensions, working at such Reynolds number s 19 i 20 . The gain in efficiency 
with respect to the FR is 20 points, whilst an order of magnitude of the maximum gain using a stator is 
typically 10 point o 19 ' 20 . 

Awaiting for more accurate local measurements of the flow angle at the exit of the CRS, a simple test 
of flow visualization with threads affixed downstream of the CRS was performed. It has been observed 
that without the RR the flow is very disorganized. When the RR is operating, at the design configuration 
(9 = 0.9 and Nfr = 2000 rpm), the flow is less turbulent, the threads are oriented with a small angle at the 
exit. This small angle seems, however to decrease when 9 is increased between 1 and 1.1. This is consistent 
with the results in section TlV Bl where it is found that the nominal operating point is observed for a value 
of 9 higher than the design value. 

The flow-rate range for which the static efficiency lays in the range 60% < r] s < 65% is: 2750 < Q < 
4150 m 3 .h _1 , that is from 76% of the nominal flow-rate up to 115% of the nominal flow-rate. One open 
question is to what extent the global performances of the CRS are affected by the axial spacing and the 
speed ratio, and whether the efficient range could be extended by varying the speed ratio. 

B. Influence of the rotation ratio 

In this paragraph, the rotation rate of FR is kept constant at Np r = 2000 rpm, and the rotation rate of RR 
is varied from to 2400 rpm. The corresponding 9 are 8 = {0 ; 0.5 ; 0.8 ; 0.85 ; 0.9 ; 0.95 ; 1 ; 1.05 ; 1.1 ; 1.15 & 1.2}. 
The axial spacing is A — 0.17. 

The overall performances of the CRS in these conditions are plotted in Fig. As expected, the more 
the rotation rate of RR increases, the more the static pressure rise of the CRS increases and the nominal 
flow-rate of the CRS increases. The maximal efficiency as a function of 9 is plotted in Fig. [5j 

For very low rotation rates of RR, i.e. for 9 — ($ in Fig. [4]) and 9 = 0.5 (< in Fig. [4j, the system is 
very inefficient: in the first case when the RR is at rest the maximum efficiency hardly reaches 35% which 
is below the maximal efficiencies of both FR and RR alone. The maximum flow-rate that can be reached is 
moreover very low in both cases compared to the discharge goal of 3600 m 3 .h _1 . 

In the range 9 € [0.8; 1.2], i.e. Nrr € [1600; 2400] rpm, the system is highly efficient. The maximum 
efficiency increases with 9 to reach a maximum value of 66.5% for 9 = 1.05 and is then quasi- const ant 
(r] s = 66.0% for 9 = 1.20). 

This is a very interesting feature of the counter-rotating system. One could imagine, simply by varying 
the RR rotation rate, to work at a constant pressure rise with an efficiency greater than 60% for a large 
flow-rate range. For instance in the present case, the system could give a constant static pressure rise of 
375 Pa with ri s > 60% for 3000 < Q < 4250 nAtr 1 with N FR = 2000 rpm, A = 0.17 and 9 6 [0.85 ; 1.2]. 

One could also imagine to work at a constant flow-rate with high static efficiency. For instance in the 
present case, the system could give a constant flow-rate of 3600 m 3 .h _1 with rj s > 60% for 290 < Ap s < 
490 Pa with N FR = 2000 rpm, A = 0.17 and 9 e [0.8 ; 1.2]. 

C. Influence of the relative axial spacing A 

Figure [S] shows the characteristics curves at the design rotation rates, i.e., Np R = 2000 rpm and 9 = 0.9. 
Regarding A e [0.17, 0.34, 0.69, 0.86], the overall performances do not change significantly and the variation 
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FIG. 4. CRS characteristics at Nfr = 2000 rpm, A = 0.17 and 8 G [0 ; 1.2] : (a) static pressure rise Ap s vs flow rate 
Q; (b) static efficiency tj s vs flow rate Q. 9 = 0, <: 9 = 0.5, >: 9 = 0.8, v= = 0.85, o: 6» = 0.9, A: 9 = 0.95, *: 
9 = 1, □: = 1.05, x: 9 = 1.1, o: = 1.15 and +: = 1.2. The blue ■ and the dashed lines stand for the design 
point of the CRS 



is in the uncertainty range. The efficiency does not vary significantly either. 

In other studies^ it was reported that the axial spacing had a more significant influence on the overall 
performances. This was noticed as well in this study. For A=2.58 and A=3.1, the global performances are 
decreased by ~ 17 Pa (5%) comparing to the other spacings. However, even for A=3.1, the CRS still shows 
good performances with high efficiency compared to the conventional fan systems. 



V. CONCLUSION 



A counter-rotating axial-flow fan has been designed according to an iterative method that is relatively 
fast. It is based on semi-empirical modelization that partly takes into account the losses, boundary layers 
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FIG. 5. Maximal static efficiency r/ s vs 9 for the CRS with N FR = 2000 rpm and A = 0.17. 



at hub and casing, and the effects of "low" Reynolds numbers (below 2 x 10 5 ). 

The overall performances at the nominal design point are slightly lower than predicted, with a static 
pressure rise 10.2% lower. The static efficiency is however remarkably high (r) s ~ 65%) and corresponds to 
a 20 points gain in efficiency with respect to the FR maximal efficiency and to a 10 points gain with respect 
to the RR. The overall measurements give first clues that allow to validate the design method. 

The counter-rotating system has a very flexible use that allows to work at constant flow-rate on a wide 
range of static pressure rises or to work at constant pressure rise on a wide range of flow-rates, with static 
efficiency bigger than 60%, simply by varying the RR rotation rate. One could thus imagine an efficient 
closcd-loop-controllcd axial-flow fan. The overall performances moreover do not significantly vary with the 
axial spacing in the range A 6 [0.17; 0.86]. However, for A = 2.58 and A = 3.1 the overall performances 
slightly decrease. 



BIBLIOGRAPHY 



1 Lesley, E., 1933. Experiments with a counter- propeller. Tech. Rep. 453, National Advisory Committee for Aeronautics. 
2 Playle, S. C, D.Korkan, K., and von Lavante, E., 1986. "A numerical method for the design and analysis of counter-rotating 

propellers". Journal of Propulsion, 2, pp. 57-63. 
3 van Gunsteren L.A., 1971. "Application of momentum theory in counterrotating propeller design" . International Shipbuilding 

Progress, 18, pp. 359-372. 

4 Chcn, B. Y.-H., and Reed, A. M., 1989. A design method and an application for contrarotating propeiiers. Tech. rep., David 
Taylor Research Center. 

5 Bechet, S., Negulescu, C, Chapin, V., and Simon, F., 2011. "Integration of cfd tools in aerodynamic design of contra-rotating 
propellers blades". In 3rd CEAS Conference (Council of European Aerospace Societies). 

6 Neise, F., and Holste, W. B., 1997. "Noise source identification in a propfan model by means of acoustical near field 
measurements". Journal of Sound and Vibration, 203, pp. 641-665. 

7 Blandcau, V. P., Joseph, P. F., and Tester, B. J., 2009. "Broadband noise prediction from rotor-wake interaction in contra- 
rotating propfans". In 15th AIAA/CEAS Aeroacoustics Conference -30th AIAA Aeroacoustics Conference, pp. AIAA 2009- 
3137. 

Peters, A., and Spakovszky, Z. S., 2012. "Rotor interaction noise in counter-rotating propfan propulsion systems". Journal 
of Turbomachinery, 134, P- 011002. 
9 Shigcmitsu, T., Furukawa, A., Watanabc, S., Okuma, K., and Fukutomi, J., 2009. "Internal flow measurement with ldv at 
design point of contra-rotating axial flow pump". Journal of Fluid Science and Technology, 4, PP- 723-734. 



11 




FIG. 6. CRS characteristics at various axial spacing: (a) static pressure rise Ap a vs flow rate Q; (b) static efficiency 
j)s vs flow rate Q. The rotation ratio of FR is Nfr = 2000rpm and 6 = 0.9. o: A = 0.17, O: A = .34, V: A = 0.69, 
x>: A — 0.86, □: A — 2.58 and X: A — 3.10. The blue ■ and the dashed lines stand for the design point of the CRS 
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